Glassy behavior of electrons near metal-insulator transitions 
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The emergence of glassy behavior of electrons is investigated for systems close to the disorder 
and/or interaction-driven metal-insulator transitions. Our results indicate that Anderson localiza- 
tion effects strongly stabilize such glassy behavior, while Mott localization tends to suppress it. We 
predict the emergence of an intermediate metallic glassy phase separating the insulator from the 
normal metal. This effect is expected to be most pronounced for sufficiently disordered systems, in 
agreement with recent experimental observations. 
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In many disordered electronic systems, electron- 
electron interactions and disorder are equally important, 
and lead to a rich variety of behaviors which remain dif- 
ficult do understand. Their competition often leads to 
the emergence of many metastable states and the result- 
ing history-dependent glassy dynamics of electrons. Such 
behavior has often been observed in sufficiently low mo- 
bility materials 0, > but more recent experiments 0, 
have provided striking and precise information of such 
non-ergodic dynamics. 

Theoretically, the possibility for glassy behavior in the 
charge sector has been anticipated a long time ago Q 
in situations where the electrons are strongly localized 
due to disorder. In the opposite limit, for well delocal- 
ized electronic wavef unctions, one expects a single well 
defined ground state and absence of glassiness. The be- 
havior in the intermediate region has proved more diffi- 
cult to understand, and at present little is known as to 
the precise role and stability of the glassy phase close to 
the metal-insulator transition (MIT) Q. On physical 
grounds, one expects the quantum fluctuations [S| asso- 
ciated with mobile electrons to suppress glassy ordering, 
but their precise effects remain to be elucidated. Note 
that even the amplitude of such quantum fluctuations 
must be a singular function of the distance to the MIT, 
since they are dynamically determined by processes that 
control the electronic mobility. 

To clarify the situation, the following basic questions 
need to be addressed: (I) Does the MIT coincide with the 
onset of glassy behavior? (2) How do different physical 
processes that can localize electrons affect the stability 
of the glass phase? In the following, we provide simple 
and physically transparent answers to both questions. 
We find that: (a) Glassy behavior generally emerges be- 
fore the electrons localize; (b) Anderson localization 
enhances the stability of the glassy phase, while Mott lo- 
calization tends to suppress it. We thus predict the 
emergence of an intermediate metallic glass phase sepa- 
rating the ordinary metal from the insulator. However, 
we expect this effect to be of relevance only for sufficiently 
strong disorder, consistent with recent experiments U 



As a simplest model where both Anderson and Mott 
routes to localization can occur, we consider an extended 
Hubbard model of disordered spin 1/2 electrons, as given 
by the Hamiltonian 

H = ^2(-tij+£i8ij)cl t(T Cj ta +U y^^iT n a+y^ VijSnjSnj. 

ija i ij 

Here, brii = m— (rij) represent local density fluctuations 
( (rii) is the site-averaged electron density), U is the on- 
site interaction, and £i are Gaussian distributed random 
site energies of variance W 2 . In order to allow for glassy 
freezing of electrons in the charge sector, we introduce 
weak inter-site density-density interactions Vij , which we 
also also choose to be Gaussian distributed random 
variables of variance V 2 / z (z is the coordination num- 
ber). We emphasize that, in contrast to previous work 
0, we shall keep the coordination number z finite, in 
order to allow for the possibility of Anderson localiza- 
tion. To investigate the emergence of glassy ordering, we 
formally average over disorder by using standard replica 
methods @, and introduce collective Q-fields to decou- 
ple the inter-site V-term 8]. As a result, the replicated 
partition function can be written in the form 

Z n = l'Y[D4Dc % DQ l de l P(e l )exp{-S}, (1) 
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where S = S el + Sq + S int , with 
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Sint = ~\v 2 Y. J J drdr'Snt(r)Qf(T-T')Snt(T'). 
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Here, Cj a (uj m ) are the replicated Grassmann 
fields Q corresponding to electrons, 5n"(r) = 
ZAc£(t)cUt)- < 4 Q (T)C(r) >], a = l,...n 
(n — ► 0) are the replica indices, and = ~/y 1 is the 
inverse lattice matrix corresponding to the interactions 
% 

Saddle-point theory. To investigate the glassy behav- 
ior on a mean-field level we formally integrate out 
the conduction electrons to produce an effective action 
S e ff[Q] for the Q- fields, the variation of which produces 
the following saddle-point conditions 

QT(r ~ r') = X(r -?*)={ KWKM) SP > (5) 
Qf(r - t') = q ab = ( 5nt(T)6nt(T>)) Sp ; (a ? b). (6) 

Physically, x( T ~ T ') is the averaged local dynamic com- 
pressibility, and q ab is related to the familiar Edwards- 
Anderson order parameter In these expressions, the 
averages are taken with respect to the saddle-point action 
for conduction electrons, as given by 

S sp = S el [clc}+S int [clc;x,q ah ] (7) 

Glass transition. In our approach, the emergence of 
many metastable states corresponding to glassy order- 
ing is identified as a replica symmetry breaking (RSB) 
instability ||. To perform such a stability analysis, we 
evaluate the S e ff[Q] at the saddle point, and then exam- 
ine its variation to infinitesimal RSB perturbations of the 
form q ab = q + Sq ab . The corresponding stability matrix 
can be expressed through appropriate correlation func- 
tion of density fluctuations, and the rest of the analysis 
is identical as in the classical case 0, 0] . The resulting 
instability criterion, which corresponds to the vanishing 
of the relevant eigenvalue of the stability matrix, takes 
the form 

\ o = l-V 2 J2\xU"n = 0)Uis=0. (8) 
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Here, the non-local static compressibilities are defined 
(for a fixed realization of disorder) as 

Xij(u n = 0) = -drii/dsj (9) 

where ni is the local expectation value of the electron 
density, and [• ■ - ]di S represents the average over disorder. 



The eigenvalue A G play the same role in our theory as 
the parameter r in a conventional Landau- Ginzburg ac- 
tion, the vanishing of which indicates an ordering insta- 
bility in the appropriate channel. Such parameters are 
typically assumed to be smooth (regular) functions of 
control parameters such as temperature or the Fermi en- 
ergy. In our case the situation is more interesting: we will 
show that the correlation function x^ = Y^,j[Xij]dis is a 
singular function, diverging at an Anderson-like metal- 
insulator transition. In this way, Anderson localization 
can be regarded as a singular perturbation in the case 
of the glassy ordering of electrons. We should emphasize 
that this unusual sensitivity to Anderson localization is 
not found in cases of more conventional transitions to 
uniform ordering. Our analysis can be easily repeated 
in such situations, and the corresponding instability cri- 
terion would instead involve the average compressibil- 
ity evaluated at the relevant ordering wave vector K , as 
given by 

x [1) (k) =X>^' r [x«K = o)W (io) 
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It is well known Q that such quantities remain finite 
(see also below) and thus non-singular at an Anderson 
transition. Similar behavior is found in the case of a su- 
perconducting instability in presence of localization |l2l | . 
In this case, the relevant pairing susceptibility (in the 
Cooper channel) was found to remain noncritical at the 
Anderson transition, opening the possibility for a direct 
superconductor-insulator transition, as seen in many ex- 
periment. 

For general values of disorder W, and interactions U 
and V, evaluating is difficult, since it has to be 

computed with respect to the action S sp [cj , c] describing 
disordered interaction electrons in finite dimensions. The 
situation is simpler both in limits of very strong and very 
weak disorder, where reliable approximations are avail- 
able. We first examine the limit of very strong disorder, 
and determine the critical value of the Fermi energy cor- 
responding to the emergence of glassy ordering at T = 0. 
To determine the transition line, to leading order in W/U 
and W/V it is sufficient to set U = V = in computing 
the required quantity ■ 

High disorder - Anderson transition. As the disor- 
der grows, the system approaches the Anderson tran- 
sition at t = t c (W) ~ W. The first hint of singu- 
lar behavior of x^ m an Anderson insulator is seen 
by examining the deeply insulating, i. e. atomic limit 
W ^S> t, where to leading order we set t = and 
obtain Xij — °~( £ i ~ We immediately find that 

X^'(K) = [5(ei — n)]dis — P(e — A*) remains finite, but 
t^( 2 ) = [5 2 (si — [i)]di8 = +°° diverges! Since we expect 
all quantities to behave in qualitatively the same fashion 
throughout the insulating phase, we anticipate x^ to 
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diverge already at the Anderson transition. Note that, 
since the instability of the glassy phase occurs already 
at %( 2 ) = V~ 2 , the glass transition must precede the lo- 
calization transition. Thus, for any finite inter-site inter- 
action V , we predict the emergence of an intermediate 
metallic glass phase separating the Fermi liquid from the 
Anderson insulator. Assuming that near the transition 

x m~A-ftiw)-B)-<* (ii) 

(A and B = t c /W are constants of order unity), from 
Eq. © we can estimate the form of the glass transition 
line, and we get 

St(W) = t G {W) — t c (W) ~ V 2/a W 1 - 2/a - 1 W -> oo (12) 

The glass transition and the Anderson transition lines 
are predicted to converge at large disorder for a < 2, 
and diverge for a > 2. Since all the known exponents 
characterizing the localization transition seem to grow 
with dimensionality, we may expect a particularly large 
metallic glass phase in large dimensions. 

In order to confirm this scenario by explicit calcula- 
tions, we compute the behavior of \^ & t the Anderson 
transition of a half-filled Bethe lattice of coordination 
z = 3. We use an essentially exact numerical approach 1 1 fill 
based on the recursive structure of the Bethe lattice [rj| . 
In this approach, local and non-local Green's functions 
on a Bethe lattice can be sampled from a large ensemble, 
and the compressibilities \ij can be then calculated by 
examining how a local charge density m is modified by 
an infinitesimal variation of the local site energy Sj on 
another site. To do this, we have taken special care in 
evaluating the local charge densities n, by numerically 
computing the required frequency summations over the 
Matsubara axis, where the numerical difficulties are min- 
imized. Using this method, we have calculated x^ as a 
function of W/t (for this lattice at half-filling Ev = 2\f2t 
), and find that it decreases exponentially [18| as the 
Anderson transition is approached. We emphasize that 
only a finite enhancement of \^ is required to trigger 
the instability to glassy ordering, which therefore occurs 
before the Anderson transition is reached. The resulting 
T = phase diagram, valid in the limit of large disorder, 
is presented in Fig. 1. Note that the glass transition 
line in this case has the form tc(W) ~ W, in agreement 
with the fact that exponential critical behavior of \ 
corresponds to a — > oo in the above general scenario. 
These results are strikingly different from those obtained 
in a theory which ignores localization 8], where ta(W) 
was found to be weakly dependent on disorder, and re- 
main finite as W — ► oo. Anderson localization effects 
thus strongly enhance the stability of the glass phase at 
sufficiently large disorder. Nevertheless, since the Fermi 
liquid to metallic glass (FMG) transition occurs at a fi- 
nite distance before the localization transition, we do not 
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FIG. 1: Phase diagram for electrons on the z = 3 Bethe 
lattice, where localization effects are treated using an exact 
numerical method. The results are valid in the limit of large 
disorder (W/U — » oo). 



expect the leading quantum critical behavior [13J at the 
FMG transition to be qualitatively modified by the lo- 
calization effects. 

Low disorder - Mott transition. In the limit of weak 
disorder W <C U, V, and interactions drive the metal- 
insulator transition. Concentrating on the model at half- 
filling, the system will undergo a Mott transition [t| as 
the hopping t is sufficiently reduced. Since for the Mott 
transition tuottiU) ~ U, near the transition W -C t, and 
to leading order we can ignore the localization effects. In 
addition, we assume that V <C U, and to leading order 
the compressibilities have to be calculated with respect to 
the action S e i of a disordered Hubbard model. The sim- 
plest formulation that can describe the effects of weak 
disorder on such a Mott transition is obtained from the 
dynamical mean- field theory (DMFT) 0. This formu- 
lation, which ignores localization effects, is obtained by 
rescaling the hopping elements t — > t/yfz and then for- 
mally taking the limit of large coordination z — > oo. To 
obtain qualitatively correct analytical results describing 
the vicinity of the disordered Mott transition at T = 0, 
we have solved the DMFT equations using a 4-boson 
method of Kotliar and Ruckenstein [J^. At weak dis- 
order, these equations can be easily solved in close form, 
and we simply report the relevant results. The critical 
value of hopping for the Mott transition is found to de- 
crease with disorder, as 

t c {W)^t°{l-4(W/U) 2 + ■■■), (13) 

where for a simple semi-circular density of states 0] 
t° = 37rl764 (in this model, the bandwidth B = At). 
Physically, the disorder tends to suppress the Mott in- 
sulating state, since it broadens the Hubbard bands and 
narrows the Mott-Hubbard gap. At sufficiently strong 
disorder W > U, the Mott insulator is suppressed even 
in the atomic limit t —* 0. The behavior of the compress- 
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FIG. 2: Global phase diagram for the disordered Hubbard 
model, as a function of the hopping element t and the dis- 
ordered strength W, both expressed in units of the on-site 
interaction U. The size of the metallic glass phase is deter- 
mined by the strength of the inter-site interaction V. 



ibilities can also be calculated near the Mott transition, 
and to leading order we find 



A 2 ) 



(1 



t c {W). 



(1 + 2%{W/U) 2 



(14) 



Therefore, as any compressibility, x is found to vanish 
in the vicinity of the Mott transition, even in presence 
of finite disorder. As a result, the tendency to glassy 
ordering is strongly suppressed at weak disorder, where 
one approaches the Mott insulating state. 

Finally, having analyzed the limits of weak and strong 
disorder, we briefly comment on what may be expected 
in the intermediate region W > U. On general grounds, 
we expect a global phase diagram as shown in Fig. 2. The 
Mott gap cannot exist for W > U, so in this region and 
for sufficiently small t (i. e. kinetic energy), one enters 
an gapless (compressible) Mott-Anderson insulator. For 
W > U, the computation of x^ requires the full solution 
of the Mott-Anderson problem. The required calcula- 
tions can and should be performed using the formulation 
of Ref. , but that difficult task is clearly beyond the 
scope of this letter. However, based on general arguments 
presented above, we expect x^ f° vanish as one ap- 
proaches the Mott insulator (W < U), but to diverge as 
one approaches the Mott-Anderson insulator (W > U). 
Near the tetracritical point M (see Fig. 2), we may ex- 
pect x (2) ~ SW- a St f3 , where 5W = W-W Mo tt(t) is the 
distance to the Mott transition line, and St — t — t c (W) 
is the distance to the Mott-Anderson line. Using this 
ansatz and Eq. Jjyl, we find the glass transition line to 
take the form 



We thus expect the intermediate metallic glass phase 
to be suppressed as the disorder is reduced, and one ap- 
proaches the Mott insulating state. Physically, glassy 
behavior of electrons corresponds to many low-lying re- 
arrangements of the charge density; such rearrangements 
are energetically unfavorable close to the (incompress- 
ible) Mott insulator, since the on-site repulsion U op- 
poses charge fluctuations. Interestingly, very recent ex- 
periments on low density electrons in silicon MOSFETs 
have revealed the existence of exactly such an interme- 
diate metallic glass phase in low mobility (highly disor- 
dered) samples 13. In contrast, in high mobility (low dis- 
order) samples |5J , no intermediate metallic glass phase is 
seen, and glassy behavior emerges only as one enters the 
insulator, consistent with our theory. Similar conclusions 
have also been reported in studies of highly disordered 
InC>2 films [3 , where the glassy slowing down of the elec- 
tron dynamics seems to be suppressed as the disorder 
is reduced and one crosses over from an Anderson-like 
to a Mott-like insulator. In addition, these experiments 
0, IU provide striking evidence of scale- invariant dynam- 
ical correlations inside the glass phase, consistent with 
the hierarchical picture of glassy dynamics, as generally 
emerging from mean-field approaches such as the one 
used in this letter. 

We thank S. Bogdanovich, S. Chakravarty, J. 
Jaroszynski, D. Popovic, and Z. Ovadiyahu for useful dis- 
cussions. This work was supported by the NSF grant 
DMR-9974311 and the National High Magnetic Field 
Laboratory. 



[1] 
[2] 

[3] 
[4] 

[5] 

[6] 



[7] 

[8] 

[9] 
[10] 



St = t G {W) - t c (W) ~ 8W 0/a ; W > W M - 



(15) 



P.A. Lee and T. V. Ramakrishnan, Rev. Mod. Phys. 57, 
287 (1985). 

M. Ben-Chorin, et al., Phys. Rev. B 48, 15025 (1993); 
Z. Ovadyahu and M. Pollak, Phys. Rev. Lett. 79, 459 
(1997). 

A. Vaknin et al, Phys. Rev. Lett. 81, 669 (1998). 

S. Bogdanovich and D. Popovic, Phys. Rev. Lett. (2002, 

in press), cond-mat/0106545 

J. Jaroszynski, Dragana Popovic, and T. M. Klapwijk, 
cond-mat/0205226 

A. L. Efros et al, J. Phys. C 8 L49 (1975); J. H. Davies, 
et al, Phys. Rev. Lett. 49, 758 (1982); Phys. Rev. B 29, 
4260 (1984); M. Pollak, Philos. Mag. B 50, 265 (1984); 
M. Pollak and M. Ortuno, Sol. Energy Mater. 8, 81 
(1982). 

N. F. Mott, Metal- Insulator Transitions (Taylor and 
Francis, London 1974). 

A. A. Pastor and V. Dobrosavljevic, Phys. Rev. Lett. 83, 
4642 (1999). 

P. W. Anderson, Phys. Rev. 109, 1498 (1958). 
Recent work 8] has shown that in presence of random 
site energies, such random interactions are generated by 
renormalization even if one starts with purely repulsive 
(e. g. Coulomb) interactions in the bare Hamiltonian. 
Similar conclusions have been obtained numerically in: 



5 



E. R. Grannan and C. C. Yu, Phys. Rev. Lett. 71, 3335 
(1993). 

[11] K. H. Fischer and J. A. Hertz, Spin Glasses (Cambridge 
University Press, Cambridge, England, 1993). 

[12] A. Kapitulnik and G. Kotliar, Phys. Rev. Lett. 54, 473- 
476 (1985). 

[13] D. Dalidovich and V. Dobrosavljevic, preprint 

|cond-ma"t /0205269 
[14] A. Georges et al. , Rev. Mod. Phys. 68, 13 (1996). 



[15] G. Kotliar and A. E. Ruckenstein, Phys. Rev. Lett. 57, 

1362-1365 (1986). 
[16] V. Dobrosavljevic and G. Kotliar, Phys. Rev. Lett. 78, 

3943 (1997). 

[17] R. Abou-Chacra, P. W. Anderson, D. J. Thouless, J. 

Phys. C 6, 1734 (1973). 
[18] A. D. Mirlin, Y. N. Fyodorov, Nucl. Phys. B 366, 507 

(1991). 



